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We present a study of B~ ^DqpK~ decays, where D^p is reconstructed in CP-even channels, 
based on a sample of 88.8 million T(45) BB decays collected with the BABAR detector at the 
PEP-II e"'"e~ storage ring. We measure the ratio of Cabibbo-suppressed to Cabibbo-favored branch- 
ing fractions S(B"-^D^piC")/B(B"^-D^p7r-) = (8.8 ± 1.6(stat) ± 0.5(syst)) x 10"^ and the CP 
asymmetry Acp = 0.07±0.17(stat)±0.06(syst). We also measure B(B--^D''A'-)/Z3(B-^_D''7r-) = 
(8.31 ± 0.35(stat) ± 0.20(syst)) x 10"^ using a sample of 61.0 million BB pairs. 

PACS numbers: 14.40.Nd, 13.25.Hw 



The recent observation of CP violation in the B meson 
system has provided a clean measurement of the angle 
/3 of the unitarity triangle. Although this measurement is 
in good agreement with the expectations of the Standard 
Model derived from other measurements of weak interac- 
tions, further measurements of CP violation in B decays 
arc needed to overconstrain the unitarity triangle and 
confirm the CKM mechanism or observe deviations from 
it. A theoretically clean measurement of the angle 7 can 
be obtained from the study of B~^D^*^°K^*'^~ decays 
by reconstructing the meson in Cabibbo-allowed CP 
eigenstates and doubly Cabibbo-suppressed decays 0]. 

In this Letter we present the measurement of the ra- 
tios of Cabibbo-suppressed to Cabibbo-favored branching 
fractions 

_ B{B-^D'{Dlp)K-) + B{B+^DO{D'cp)K+) 
{CP) B[B-^D°{Dlp)n-) + B{B+^D°{DIp)tt+) 

(1) 

with D° reconstructed in Cabibbo-allowed or CP-even 
(D^p) channels. The direct CP asymmetry 

^ BjB-^D^cpK-) - B{B+^D°cpK+) 

BiB-^D°cpK-) + BiB+^D°cpK+) 

is also measured. The measurement of R uses a sample 
of 61.0 million T(4S') decays in BB pairs collected with 
the BABAR detector at the PEP-II asymmetric-energy B 
factory. The analysis of B^^D'(jpK^ decays uses a sam- 
ple of 88.8 million BB pairs. Since the BABAR detector 
is described in detail elsewhere only the components 
that are crucial to this analysis are summarized here. 
Charged-particle tracking is provided by a five-layer sil- 
icon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH). For charged-particle identification, ionization en- 
ergy loss in the DCH and SVT, and Cherenkov radia- 
tion detected in a ring-imaging device (DIRC) are used. 
Photons are identified by the electromagnetic calorimeter 
(EMC), which comprises 6580 thallium-doped Csl crys- 
tals. These systems are mounted inside a 1.5-T solenoidal 
superconducting magnet. We use the GEANT ^ soft- 
ware to simulate interactions of particles traversing the 
detector, taking into account the varying accelerator and 
detector conditions. 

We reconstruct B^^D'^h^ decays, where the prompt 
track h~ is a kaon or a pion (reference to the charge- 
conjugate state is implied here and throughout the text 
unless otherwise stated). Candidates for I?° are re- 
constructed in the non-CP flavor eigenstates K^tt^, 



K TT+TT+TT , K tt+tt" (nou-CP modes) and in the CP- 
even eigenstates 7r^7r+ and K^K^ {CP modes). 

To reduce the combinatorial background, only charged 
tracks with momenta greater than 150 MeV/c are used 
to reconstruct D'^^if^Tr+Tr+Tr" and D'^-^K~7r~^Tr'^; the 
prompt particle h is required to have momentum greater 
than 1.4 GeV/c. Particle ID information from the drift 
chamber and, when available, from the DIRC must be 
consistent with the kaon hypothesis for the K meson can- 
didate in all D'^ modes and with the pion hypothesis for 
the TT^ meson candidates in the P)°— >7r~7r+ mode. For 
the prompt track to be identified as a pion or a kaon, we 
require that its Cherenkov angle be reconstructed with 
at least five photons. We reject a candidate track if its 
Cherenkov angle is consistent with that of a proton or if 
it is identified as an electron by the DCH and the EMC. 

Photon candidates are required to have energies grater 
than 70 MeV. Photon pairs with invariant mass within 
the range 124-144 MeV/c^ and total energy greater than 
200 MeV are considered tt" candidates. To improve the 
momentum resolution, the 7r° candidates are kinemati- 
cally fit with their mass constrained to the nominal tt*' 
mass 0. 

The invariant mass of a I?" candidate, A/(P'°), must 
be within 3a of the mean fitted mass for the channels 
K^TT^, /^"Tr+Tr+Tr" and K~K~^, and within 2a- for the 
K^tt^tt'^ channel. Candidates for D^^tt^tt^ are se- 
lected in the range 1.80 < M{D°) < 1.93 GeV/c^ and 
the invariant mass of the (/i~7r"'") system, where tt"*" is the 
pion from and h~ is the prompt track taken with the 
kaon mass hypothesis, must be greater than 1.9 GeV/c^ 
to reject background from B~^D^[-^ if~7r+]7r~ and 
B^^K*[-^ X^tt+Jtt^ decays. For aU the D° decay 
channels except the tt^tt"*" mode a kinematical fit to the 
nominal I?" mass |5i is applied. The Z?°— >7r~7r"'" selection 
differs because of its particular background, as described 
later. 

We reconstruct B meson candidates by combining a 
Z?" candidate with a track h. For the non-CP modes, 
the charge of the track h must match that of the kaon 
from the meson decay. We select B meson candi- 
dates by using the beam-energy-substituted mass toes = 
^{E'f- jl + pj • ipsY /Ef - p| and the energy difference 
/S.E — Eg — E* /2, where the subscripts i and B refer to 
the initial e~^e~ system and the B candidate respectively, 
and the asterisk denotes the CM frame. The toes distri- 
butions for B~^D'^h~ signals are Gaussian distributions 
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centered at the B mass with a resolution of 2.6MeV/c^, 
which does not depend on the decay mode or on the 
nature of the prompt track. In contrast, the AE dis- 
tributions depend on the mass assigned to the prompt 
track and on the D'^ momentum resolution. We evaluate 
AE with the kaon mass hypothesis so that the distri- 
butions are centered near zero for B~—>D'^K~ events 
and shifted by approximately 50MeV for B~ ^D'^tt~ 
events. The AE resolution is about 20MeV for the 
>7r~7r+ mode, and typically 17MeV for the other 
D° decay modes. We select B mesons in the range 
5.2 < rriEs < 5.3 GeV/c^ with the exception of the 
D''^7r~7r"'"mode, for which toes is required to be within 
3(7 of the mean value. All B candidates are selected in 
the range -0.10 < AE < 0.13 GeV. For events with mul- 
tiple _B~— >D''/i~candidates, the best candidate is chosen 
based on the values of M(Z?°) and toes; this happens 
in fewer than 1% of the selected events for two-body D'~' 
decays and in « 4% of the events for the other D° decays. 

To reduce backgrounds from continuum production of 
light quarks, we make use of two quantities that exploit 
the different topologies of e+e" —^qq{q = u, d, s, c) and 
BB events. The first quantity is the normalized sec- 
ond Fox- Wolfram moment i?2 = H2/HQ, where Hi 
is the ^-order Fox- Wolfram moment of all the charged 
tracks and neutral clusters in the event. Only events 
with i?2 < 0.5 are selected. The second quantity is the 
angle 9t between the thrust axes of the B candidate 
and of the remaining charged tracks and neutral clus- 
ters, evaluated in the CM. We require |cos0t| < 0.9 
for the D^^K~-K'^ mode, and \ cos6t\ < 0.7 for the 
D^^K-TT+TT+TT- and D°^K-Tr+TT° modes. For the 
D^^K^ and £'"^7r^7r+ modes an additional quan- 
tity is used to suppress further the continuum back- 
ground: the angle Ooh between the direction of one of 
the decay products of the 13° and the direction of flight 
of the B, in the rest frame. The quantities cos Ot 
and cos9]jh are uncorrelated for the signal but not for 
the continuum background. This correlation is exploited 
to make a more efficient cut in the cos 9t — cosOuh plane. 

The total reconstruction efficiencies, based on simu- 
lated signal events, are 42%(ii'"7r+), lA%{K^Tr+7r+Tr^), 
8%iK-n+7r°), 3A%{K-K+) and 36%(7r-7r+). 

The main contributions to the BB background for the 
non-CP modes come from the processes B^D*h [h = 
TT, if), B~^D^p~ and mis-reconstructed B~^D^h~'. 
For D^p decays, the backgrounds B~^K~K~^K~ and 
B^ ^K^TT^TT" 7] must also be considered, since they 
have the same AE and toes distribution as the D'^K^ 
signal. The resonant component of these decays is negli- 
gible after the selection requirements for the CP modes. 

For each decay mode an extended unbinned 
maximum-likelihood fit to the selected data events deter- 
mines the signal and background yields rii [i — 1 to M, 
where M is the total number of signal and background 
channels). Two kinds of signal events, B~ ^D^-k~ and 



B~^D^K~ , are considered, while the number of back- 
ground sources depends on the channel. For non-CP 
modes we consider four kinds of backgrounds: candidates 
selected either from continuum or from BB events, in 
which the prompt track is either a pion or a kaon. In the 
case of D^^K^ we consider two kinds of background 
depending on the nature of the prompt track. Finally, in 
the case of D'^ ^■k~ tt^ we consider four contributions: 
the B^— ^if ^TT+TT^ and B^ tt^tt+tt^ decays and two 
kinds of generic background depending on the nature of 
the prompt track. 

The input variables to the fit for the non-CP and the 
D^^K~ K'^ modes are toes, AE, and a particle identi- 
fication probability for the prompt track based on the 
Cherenkov angle 9c, the momentum p and the polar 
angle 9 of the track. For the D^^tt'tt^ mode, toes 
is replaced by M(D'^). This allows us to separate the 
B^ from the non-resonant B^ ^K^tt^tt^ con- 
tributions since the (7r^7r+)-invariant-mass distribution 
peaks at the D'^ mass for signal while it is featureless 
for background. The extended likelihood function C is 
defined as 



M 



N 



exp 



-E- n 



M 



(3) 



where N is the total number of observed events. The M 
functions Vi (xj ; di ) are the probability density functions 
(PDFs) for the variables Xj , given the set of parameters 
di. They are evaluated as a product Vi — 'Pi{AE,x) x 
'Pi{9c) where x = toes or M{D^) depending on the 
channel. 

The Gaussian shape of the toes PDF for signal 
events is determined from a pure sample of _B~— >D*'7r~, 
D'^^K^tt^ decays selected from on-resonance data. The 
AE distribution for B^^D'^K^ signal events is param- 
eterized with a Gaussian distribution whose parameters 
are determined from a pure sample of B" ^D^Tr~ events 
selected after the pion mass is assigned to the prompt 
track. The displaced AE distribution for B^— s-D^tt" is 
parameterized with a sum of two Gaussian distributions. 
The shape of M{D°) is also described by a Gaussian dis- 
tribution whose parameters are determined from data. 

The parameters for the AE and toes distributions 
for the continuum background of the non-CP and 
D^^K~ K'^ modes are determined from off-resonance 
data. The background shape in AE is parameterized 
with a linear function, while that of the toes is parame- 
terized with an ARGUS threshold function [8j /(toes) oc 
toes\/1 - 2/^exp[-^(l - y"^)], where y = toes/toq and 
Too is the mean CM energy of the beams. The corre- 
lation between toes and AE for the generic BB back- 
ground is taken into account with a two-dimensional PDF 
determined from simulated events through a method 
based on the Kernel Estimation technique. For 
the D^-^K^ mode the contribution from the non- 
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resonant B~^K~K~K^ decays is estimated and 
added to the continuum and the generic BB background. 

The AE and M{D^) distributions for the continuum 
and generic BB background of the D^^tt^tt'^ mode 
are determined from off-resonance data and simulated 
events, respectively. The AE distribution is described by 
a linear function while M{D'^) is parameterized with the 
sum of a linear function (combinatorial background) and 
a Gaussian distribution (real D°^7r^7r+). The M{D^) 
PDFs of the non-resonant B~^h~Tr~7T~^ decays are de- 
scribed by linear functions, while the AE distributions 
are parameterized with one or two Gaussian distribu- 
tions, as for the B~ ^D'^h~ signals. 

Finally, the parameterization of the particle identi- 
fication PDF is performed by fitting with a Gaussian 
distribution the background-subtracted distribution of 
the difference between the reconstructed and expected 
Cherenkov angles of the kaons and pions from decays, 
in a pure D*'^ — > Z^^tt"*" {D^—^K~ti^) control sample. 

The results of the fit are summarized in Table ^ Fig- 
ure n shows the distributions of AE for the combined 
non-CP and CP modes after enhancing the B — » D^K 
purity by requiring that the prompt track be consistent 
with the kaon hypothesis and that |mES — (wes)| < 3ct 
(|M(7r-7r+) - (M(7r-7r+))| < 3cr for D°^tt-tt+). The 
projection of a likelihood fit, modified to take into ac- 
count the tighter selection criteria, is overlaid in the Fig- 
ure. 



TABLE I: Results from the maximum-likelihood fit. For the 
D^^~*K-K+ and — >7r TT modes we quote the results for 
the fits performed on the whole sample and on the B'^ and 
B~ subsamples. 





mode 


N{B D"n) 


N{B D"K) 


N{r{4S)) 


K 


-7V + 


4440 ± 69 


360 ± 21 


61.0 X 10" 


K~ 




2914 ± 56 


242 ± 18 


61.0 X 10"^ 


K' 




2650 ± 56 


208 ± 18 


61.0 X 10*^ 


K- 


-K+ 


565 ± 25 


44.3 ±9.0 


88.8 X 10" 


K- 


-K+ [B+] 


286 ± 18 


16.7 ±5.8 


88.8 X 


K- 


-K+ [B-] 


280 ± 18 


27.8 ±6.8 


88.8 X 10"^ 


TX~ 


7r+ 


195 ± 17 


24.2 ± 7.2 


88.8 X 10" 


■K~ 


7r+ [B+] 


99 ± 12 




88.8 X lO'' 


7r~ 


n+ [B-] 


96 ± 12 




88.8 X lO'' 



The ratios R and Rcp are computed by scaling the 
ratios of the numbers of B^ ~*D^ K~ and B~ —*D^tt^ 
mesons by correction factors that account for small 
differences in the efficiency between B~ and 
_B~^_D'^ 7r~ selection, estimated with simulated signal 
samples. The results are listed in Table Ull 

The direct CP asymmetry Acp for the B^ 
decays is calculated from the measured yields of positive 
and negative charged meson decays reported in Table 
We measure Acp = 0.07 ± 0.17(stat) ± 0.06(syst). 




0.05 0.1 

AE (GeV) 

FIG. 1: Distributions of AE for events enhanced in B — > D'^K 
signal. Top: D" K~ , K~ n^Tv'^iv~ , K~ n^n'^; bottom: 
, TT TT . Solid curves represent projections of 
the maximum likelihood fit; dashed, dashed-dotted and dot- 
ted curves represent the B D^K, B D"iy and back- 
ground contributions. 



Systematic uncertainties in the ratios R, Rcp and in 
the CP asymmetry Acp arise primarily from uncertain- 
ties in signal yields due to imperfect knowledge of the 
PDF shapes. The parameters of the analytical PDFs are 
varied by ilu and the difference in the signal yields is 
taken as a systematic uncertainty. When a BB PDF is 
parameterized through the kernel estimation, we repeat 
the fit using several statistically independent simulated 
BB samples to define the PDF. The width of the distri- 
bution of the difference between the new yields and the 
original yield is taken as the systematic uncertainty. 

The uncertainties in the branching fractions of the 
channels contributing to the BB background have been 
taken into account. The correlations between the differ- 
ent sources of systematic errors, when non- negligible, are 
considered. An upper limit on intrinsic detector charge 
bias due to acceptance, tracking, and particle identifi- 
cation efficiency has been obtained from the measured 
asymmetries in the processes B~ ^D^[^K~Tr^]h~ and 
B~ ^DQp7r~ , where CP violation is expected to be neg- 
ligible. This limit (0.04) has been added in quadrature to 
the total systematic uncertainty on the CP asymmetry. 

In conclusion, we have reconstructed B~—^D'^K~ de- 
cays with D° mesons decaying to non-CP and CP-even 
eigenstates. The ratios R^cp) of the branching frac- 
tions B{B-^D'^^p^K-) and B{B-^Df(^p^Tr-) and the 
direct CP asymmetry Acp have been measured. The 
measured ratio R is consistent with Standard Model 
expectation (« 7.5%) assuming factorization ^ld\. In 
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TABLE II: Measured ratios R and Rcp for different D decay 
modes. The first error is statistical, the second is systematic. 



B -^D"h decay mode 


B{B^DK)/B(B^D-k) (%) 




8.4 ±0.5 ±0.2 




8.7 ±0.7 ±0.2 




7.7 ±0.7 ±0.3 


weighted mean 


8.31 ±0.35 ±0.20 




8.0 ± 1.7 ±0.6 




I2.9±4.0ti-5 


weighted mean 


8.8 ± 1.6 ±0.5 



the Standard Model Rcp/R = 1 ± ± 2r cos (5 cos 7 
and Acp — 2r sin 5 sin 7/(1 ± ± 2r cos (5 cos 7), where 
r w 0.1 — 0.2 is the magnitude of the ratio of the ampli- 
tudes for the processes B~ — > D^K^ and D^K^ , 
and 5 is the (unknown) relative strong phase between 
these two amplitudes 0- The measured values of R and 
Rcp are equal within errors, and Acp is consistent with 
zero. These results, together with the ones obtained by 
CLEO and Belle represent a first step towards the 
measurement of the angle 7 and of direct CP violation in 
the B system using the B^ ^D^K^ decays. 
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